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[Plate  9] 


Tho  air  velocity  has  been  measured  in  blast  waves  produced  by  the  detonation  of  trinitro¬ 
toluene  (t.n.t.)  charges  varying  in  mass  from  30  to  200000  lb.  The  technique  consists 
essentially  of  using  a  high-speed  camera  to  record  the  displacement  of  smoke  trails  formed 
close  to  the  charge  just  before  detonation.  The  initial  decay  of  velocity  behind  the  shock  agrees 
well  with  theoretical  predictions,  such  as  those  of  Brode  (1959).  but  at  later  times  them 
is  an  extended  outward  flow,  which,  it  is  postulated,  is  caused  by  the  ‘after-burning’  of  the 
detonation  products  in  the  presence  of  atmospheric  oxygen.  It  has  been  shown  that  this 
phenomenon  does  not  occur  in  the  case  of  the  detonation  of  an  explosive  with  a  high  oxygen 
balance,  or  for  a  nuclear  detonation.  The  velocity  decay  within  a  t.n.t.  blast  wave  noay  be 
described  by  the  equation 

V  »  F’,(l-;»)oxp(-a<)-|-oln(l-l-^) 

fltted  to  the  data  by  an  iterative  least  squares  procedure.  It  has  been  demonstrated  that 
particle  trajectories,  determined  by  the  smoke  tracer  technique,  may  be  used  to  calculate  the 
variations  of  density,  pressure  and  temperature  within  the  wave,  without  reference  to 
other  measurements. 


1.  Introduction 

Sir  Geoffrey  Taylor  (1950)  calculated  the  basic  characteristics  of  an  intense  point 
source  explosion  of  negligible  mass,  and  was  able  to  describe  the  resulting  blast 
wave  to  the  time  at  which  the  peak  hydrostatic  pressure  had  fallen  to  about  10  atm. 
At  lower  pressures  his  analysis  ceased  to  be  accurate  because  the  similarity  assump¬ 
tions  that  he  had  used  to  describe  the  blast  wave  were  valid  only  when  the  peak 
pressure  was  large  compared  to  the  ambient  atmospheric  pressure.  Von  Neumann 
&  Bethe  ( 1 958)  made  similar  calculations  and,  in  addition,  gave  a  mainly  qualitative 
description  of  the  point  source  blast  wave  in  intermediate  and  low  pressure  regions 
by  assuming  in  the  first  case  that  y-l  was  small,  and  for  the  low-pressure  region 
that  an  approximate  acoustical  theory  could  be  applied.  This  theoretical  work  has 
been  further  extended  by  Sakurai  (1953,  1954)  and  Oriiima  (i960).  The  advent  of 
high-speed  digital  computers  made  it  possible  to  carry  out  complete  numerical 
solutions  for  the  spherically  symmetrical  explosion.  Results  of  such  analyses  were 
reported  by  Broyles  ( i960)  and  Brode  (1955)  for  a  point  source,  and  Brode  has  made 
similar  calculations  fw  the  explosion  of  high-pressure  spheres  (Brode  1956)  and 
spherical  trinitrotoluene  (t.n.t.)  chargee  (Brode  1959).  These  calculations  are  par¬ 
ticularly  valuable  beoause  they  give  not  only  the  shook  fifont  conditkms  but  also 
the  predicted  variation  of  piesfure,  density,  temperature  and  particle  velocity 
throughout  the  r^;ion  firom  the  charge  centre  to  the  shock  front. 

In  19C6,  a  programme  started  at  Suffield  Experimental  Station,  Canada, 
to  study!  the  physios  of  blast  waves,  using  mainly  charges  of  t.n.t.  ranging  in  site 
from  8  to  lO^lb.  tJntQlHat  time  the  blast  parameter  most  commonly  measured 
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haU  lH>en  the  hydrostatic  oveqtressun*.  The  peak  value  could  Iw  dctcnninod  from 
the  shock  velocity,  and  the  decay  with  time  hy  means  of  pieio-electric  transducers. 
Many  of  the  effects  associated  with  blast  waves  arise  from  the  dynamic  pressure  cf 
the  air  flow  after  the  passap^o  of  the  shock  front,  ami  for  this  rcasot\  studies  wore 
started  to  determine  the  air  density  and  velocity  in  this  region.  A  gauge  which  has 
been  developed  to  determine  the  density  has  been  described  by  Dewey  ft  Anson 
( •9^3)'  techniques  and  results  of  the  velocity  measurements  are  the  subject 

of  this  paper. 

2.  ExPEKIMRNTAL  TKCHNIQUBS 

Initial  studies  of  blast  waves  at  Suffield  Experimental  Station  were  concentrated 
on  the  region  in  which  the  peak  side-on  overpressure  varies  from  about  7-5  to 
O'fl  atm.  This  is  the  region  of  most  interest  to  those  concerned  with  the  rcs|M)nse  of 
structures  to  blast  waves.  At  peak  overpressures  of  7*5  and  0-5  atm  the  Ranldno  - 
Hugoniot  equation  yields  peak  particle  velocities  of  approximately  Mach  2  and 
0-3  respectively,  relative  to  the  speed  of  sound  ahead  of  the  shook,  and  these  values 
might  be  expected  to  decay  to  sero  in  times  of  about  10  to  200  ms  for  the  range  of 
explosive  yields  considered.  The  violent  accelerations  involved  have  precluded 
the  use  of  any  normal  type  of  anemometer,  although  StrOm  (196a)  has  described 
a  corona  discharge  absolute  anemometer  which  might  be  used  for  this  purpose. 
The  technique  developed  at  HufHeld  uses  high-speed  cameras  to  record  the  move- 
mcntofvisihlc  low  inertia  tracers.  In  conjunction  with  this  work,  Muirhend.Lecuyer 
ft  McCaliiim  (1959)  and  Muirhead  ft  Lccuycr  (1959).  measured  the  movement  of 
various  types  of  smoko  streams  in  shock  tube  flows,  and  within  the  limits  of  experi¬ 
mental  error  showed  that  they  acted  as  perfect  tracers  of  the  air  flows.  This  tech¬ 
nique  was  later  extended  to  study  the  air  flows  associated  with  the  interaction  of 
shock  waves  with  a  variety  of  small-scale  structures  (Muirhead  ft  McCallum  1959). 

The  method,  as  developed  for  use  in  the  field  with  large  yield  explosions,  con¬ 
sisted  of  forming  smoke  trails  in  a  vertical  plane  containing  the  charge  centre  and 
parallel  to  the  object  plane  of  a  high-speed  cin<  camera.  The  displacement  with  time 
of  the  trails  was  measured  from  the  resulting  films;  the  air  velocity  was  then  calcu¬ 
lated,  and  subsequently  transformed  from  Lagrangian  to  Eulerian  co-ordinates  to 
give  the  velocity-time  variation  at  specific  distances  from  the  explosion.  A  typical 
set  of  smoke  trails  is  shown  in  figure  1,  plate  0,  which  is  a  photograph  taken  1 10  ms 
after  the  detonation  of  a  2000001b.  surface  burst  hemispherical  charge  of  t.n.t. 

The  trails  were  formed  by  a  ehloro-sulphonic  acid  mixture  released  from  spinning 
mortar  shells  consisting  essentially  of  steel  tubes  scaled  at  both  ends,  as  shown  dia- 
grammatically  in  figure  2.  Two  holes  in  the  walls  of  the  tube,  close  to  the  leading  end 
are  sealed  internally  by  a  metal  plug  silver  soldered  into  position.  A  7  g  hammer  is 
held  by  a  brass  shear  wire  in  a  hole  drilled  into  the  threaded  plug  which  seals  the 
shell.  At  the  lower  end  there  is  an  impulse  turbine.  The  shell  is  fired  by  means  of  a 
25  g  gunpowder  charge  contained  in  a  cartridge,  and  the  resulting  explosion  gives 
the  shell  an  initial  velocity  of  about  700  ft/s  and  an  angular  velocity  of  about 
400rev/s.  The  hammer  at  the  upper  end  is  released  by  the  impulsive  motion  of  the 
shell,  and  strips  off  the  metal  plup.  The  deceleration  by  air  resistance,  and  the 
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oentrifuf;al  forces  of  the  spinning  shell  ejects  the  30  ml.  of  acid  mixture  through 
the  unplugged  holes.  Tae  acid  rapidly  absorbs  atmospheric  water  vapour  to  form 
an  acid  droplet  trail  about  2  ft.  wide.  The  mortars  are  fired  a  few  tenths  of  a  second 
before  the  charge  is  detonated  to  produce  trails  of  the  desired  height  immediately 
before  the  arrival  of  the  shook  front. 

In  order  to  determine  accurately  distances  in  the  plane  of  the  smoke  trails, 
marker  boards  were  placed  at  intervals  along  the  mortar  line.  These  were  made 
from  steel  plates  supported  by  girders  set  in  concrete,  so  that  they  were  not  dis¬ 
placed  by  the  blast  wave,  'fhe  surface  of  each  plate  was  painted  with  a  black  and 
white  pattern  and  the  distance  and  elevation  of  the  pattern  centre  was  surveyed 
relative  to  the  charge  oenti'e. 


Ktoviuc  2.  Diagram  of  amoke  mortar  sbell. 

In  the  majority  of  experiments  two  16  mm  cameras  were  used  to  record  the 
smoke  movement;  a  Fastax  WFST  operating  at  between  2000  and  4000  frames/s, 
and  a  high-speed  Kodak  operating  in  the  range  from  1000  to  3000  fritmes/s.  Neon 
bulbs  within  the  cameras  marked  the  edges  of  the  films,  using  a  pickup  from  the 
detonation  pulse  to  obtain  a  sero  signal,  and  the  output  from  a  crystal  oscillator  to 
give  riming  marks  at  1  mr  2*6  ms  intervals. 

AU  the  experiments  described  here  were  carried  out  in  bright  sunny  conditions, 
and  it  was  found  that  Uie  contrast  of  the  smoke  trails  against  a  clear  or  hasy  sky 
was  so  slight  that  the  trails  could  be  studied  only  during  a  short  period  in  which 
they  wore  illuminated  by  the  fireball.  In  order  to  improve  the  contrast,  riic  trails 
were  photographed  against  a  black  smoke  background  fmmed  by  burning  petrol, 
gelled  by  adding  2%  Octoic  acid  and  4%  Octal  powder. 

3.  Ak ALYTICAL  PBOCEDUaia 

Ihe  individual  frames  of  the  films  were  projected  at  a  magnification  of  i^iproxi- 
matelly  20.  The  diaige  centre  was  not  normally  b  the  field  of  view,  but  its  porition 
in  the  projection  plane  could  be  determined  from  those  of  the  marker  bouds.  A 
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grid  of  radial  lines  was  drawn  on  the  projection  plane  through  this  point.  Each 
frame  was  projected  onto  the  grid  and  alined  correctly  by  positioning  the  images  of 
the  markers.  A  typical  grid  and  frame  is  shown  as  figure  3. 

It  was  assumed  that  the  air  movement  was  radial  so  that  the  intersection  of  a 
smoke  trail  and  a  grid  line  defined  a  smoke  element,  the  motion  of  which  could  be 
followed  from  frame  to  frame.  The  assumption  of  radial  movement  has  been  con¬ 
firmed  in  several  experiments  in  which  intersecting  patterns  of  smoke  trails  in  the 
same  plane  were  used,  so  that  the  movement  of  each  intersection  could  be  ob  served. 


^camera  Held  of  virw 

-  -- 

Fiuuaa  S.  Smoke  trail  pattern  projected  on  to  radial  grid.  The  intereection  of  a  trail  and  a 
grid  line  defines  a  emoke  element  the  diaplaoement  of  which  from  the  charge  centre  may 
be  measured.  The  ceineca  field  of  view  is  that  used  for  a  10000  lb.  t.n.t.  experiment  with 
trails  formed  at  a  nominal  elevation  of  OO'".  The  centre  marker  board  was  180  ft.  from 
the  charge  centre. 


In  all  oases  of  ground  burst  hembpherioal  chargee  the  movement  was  found  to  be 
radial  over  the  region  of  dispiaoement  of  a  single  mtenection.  A  comjfiete  pattern 
of  intttseoting  traUs  was  not  normally  used  in  these  experiments  since  it  entailed  a 
much  wider  grid  spacing  and  a  corresponding  reduction  in  the  quantity  of  data 
obtained.  The  technique  has  been  used,  however,  in  cases  whMe  the  flow  was  not 
expected  to  be  radial,  such  as  for  venting  undmgiottnd  exidosions. 

The  horixontal  distaitoe  of  each  element  from  the  ehai^  centre  was  measured 
with  the  aid  of  a  film  reader  consisting  essentially  of  a  potentiometer  giving  a  digital 
output  on  pundted  paper  tape.  References  labdling  the  trail,  grid  line  and  the  thne 
after  charge  detonatkm  wwe  also  coded.  The  diaplaoemmt-tlme  measurements  of 
each  ^mient  were  then  sorted  from  the  total  data  by  means  of  a  digital  ocuDputer. 
At  this  stage  the  horixontal  component  of  an  demmt  displacement  was  corrected 
for  the  effects  oi  wind,  and  resdi^  to  give  its  radhd  displacement  from  the  diatge 
centre.  In  eadi  experiment  the  wind  speed  and  direction  were  monitored  continu¬ 
ously  sA  a  limited  number  positions  and  he^ts.  Owing  to  spatial  variatkwf, 
however,  it  has  been  found  frur  more  reliable  to  determine  the  qpedfic  wind 
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component  in  the  plane  of  the  trails  by  neasuring  their  drift  before  the  arrival 
of  the  shock.  No  experiments  were  carried  out  with  winds  of  greater  than  about 
lOmiles/h. 

A  fourth-degree  polynomial  was  then  fitted  to  the  displacement-time  data  for 
each  smoke  element.  The  fit  was  made  only  up  to  the  arrival  of  the  secondary  shock, 
at  which  point  there  is  a  discontinuity  of  the  first  derivative.  Figure  4  shows  a 
typical  set  of  data  and  the  resulting  fitted  curve. 


Fiuurb  4.  Variation  of  particfe  diaplaoement  with  time.  Points  are  the 
obaen-ixi  valuta.  Line  is  the  fit  to  a  fourth -deiCCtee  ptrfynomial. 

The  standard  error  of  such  fits  varied  from  0*08  ft.  for  60  lb.  charges  to  O'Sft.  for 
a  100-ton  charge,  and  these  figures  are  probably  the  best  indication  of  Uie  accuracy 
with  which  a  trail  centre  could  be  defined.  FVom  this  point  of  the  analysis  the  para¬ 
meters  were  scaled  as  follows.  Radial  distances,  r,  were  reduced  to  r/<S  and  times, 
f,  to  etIS,  where  S  m  WP^P)-,  c  was  the  ambient  velocity  of  sound,  W  the  charge 
mass  in  pounds,  P  the  atmospheric  pressure,  and  P^  a  standard  pressure,  taken  for 
the  altitude  of  SuflSeld  to  be  13*67  Ib./in.'  or  as  14*7  at  sea  level.  The  validity  of  this 
scaling  procedure  is  discussed  in  |6.  At  0*1  ft.  intervals  of  sealed  distance,  within 
the  range  covmd  by  the  movement  of  a  specific  snuAe  element,  the  fitted  displace¬ 
ment-time  equation  was  solved  and  differentiated  to  give  a  series  of  distance-time* 
veloeity  co-ordinates.  The  resulting  co-ordinates  for  all  the  smoke  elements  were 
sorted  so  as  to  gather  the  velority-time  palm  for  each  of  the  specified  distanoes. 
A  typical  set  of  sudi  data  obtained  in  an  experiment  with  a  SOOOOOIb.  surface 
burst  charge  is  plotted  in  figure  6  where  the  given  times  are  measured  firom  the  arrival 
of  the  ahorit  at  that  distance. 

Die  error  (d'akast  squares  fit  is  normally  greatest  at  the  limits  of  the  data  and  so 
the  derivatives  at  the  extremes  (d’the  fitted  dis|daceiiieot-time  curves  were  lejected. 
Values  of  the  particle  velocity  immediately  behind  the  abodt  were  not  rddained 
ihetefure.  In  order  to  derive  this  peak  value  it  was  necessary  to  extrapolate  the 
velocity-time  oo-ordinatee  to  the  time  of  sbodt  airivai  and  thfai  oonld  be  done  most 
oonsistetitly  snd  conveniently  by  fitting  sn  equation  to  these  values.  This  had  the 
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added  advantage  that  the  resulting  equation  could  be  used  to  describe  the  shape  of 
the  velocity  decay  in  a  ooncise  form. 

The  equation  which  has  been  used  most  commonly  to  describe  the  decay  of  blast 
waves  is  a  modified  form  of  an  equation  discussed  by  Friedlander  ( 1946),  namely 

P-P.fl-f/rle-', 

where  P  is  the  blast  parameter  under  discussion,  P,  its  value  immediately  behind  the 
shock,  t  the  time  measured  from  the  shook  arrival,  a  a  decay  constant,  and  r  the 
positive  duration  at  which  time  the  parameter  first  returns  to  its  ambient  value. 


Fiui'RK  5.  VarisUon  of  >'viooit,v  with  time  at  a  Mail'd  dwtaiwf*  of  T-S  fi.  S. 

Pointa  ( • )  arv  the  obwyved  v^uca,  aoUd  line  is  the  At  to 

broken  line  ia  the  variation  predicted  from  BroileV  calctilatiooa. 

It  can  be  shown  that  this  equation  may  be  used  to  deacribe  the  di‘cay  of  aonw'  blast 
parameteri  from  t.n.t.  exfdaskms  within  limited  regions,  but  when  attempts  were 
made  to  fit  the  velocity  data,  it  became  clear  that  tltere  was  a  significant  departure 
from  the  Friedlander  form.  A  good  tit  could  mnially  bo  obtained  in  the  region  ctose 
to  the  shock  and  this  gave  a  value  of  peak  velocity  which  agreed  weli  with  the  value 
derived  from  the  shock  vdodty.  Fits  could  niA  be  obtained  at  later  times.  Any 
rsasonable  frNm  of  the  equation  significantly  undermtimated  the  velodtiea  in  this 
period,  and  it  became  dear  that  an  extended  form  of  the  equation  would  be  required 
to  deacribe  completely  the  velocity  decay.  Home  of  the  equations  investigated 

F-a  +  ^e--, 

V  -  i;fae  •'  +  (1  -ole-’«l(l  -fit). 

F-  i;(l-^)e 


were 


and 
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where  a.  A,  y  and  n  were  additional  fitted  parameters.  Consideration  of  the 
difference  between  the  observed  and  the  predicted  flow,  which  is  discussed  f)elow 
(see  §6),  indicated  an  additional  Ic^arithmic  term,  and  the  equations 

K-i;(l-y?f)e-“<+aln(l  +  f), 

V  =  f'(l->?f)e-*'  +  aln(l+yf), 

and  F  =  ^'(l-/?<)e-“'  +  aln(l+/f^) 

were  r’’  fitted,  each  giving  a  progressively  better  fit,  as  judged  by  the  standard 
error.  The  excellent  way  in  ”/hich  this  last  equation  fits  the  data  in  figure  6  is 
typical  of  that  throughout  tl  o  rej^ion  studied,  and  for  all  charge  masses  that  have 
been  used. 


4.  Cl  ?tVE-FITTINQ  METHODS 

Equations  such  as  those  discussed  above  cannot  be  fitted  to  experimental  data 
by  usual  least  squares  methods  since  the  normal  equations  contain  the  relevant 
coefficients  in  non-linear  form.  An  iterative  technique  was  therefore  developed, 
based  on  a  method  of  Gauss  and  described,  among  others,  by  Whittaker  &  Robinson 
(1924),  and  Moore  &  Zeigler  (i960).  Consider  a  function  y  =f{a,b,c,x)  vhere  x 
is  the  independent  variable  and  a,  b  and  c  are  the  parameters  for  which  it  is  required 
to  obtain  a  best  fit.  If  =  £[^— /(c,6,c)]*,  where  p  is  the  observed  dependent 
variable,  it  is  required  to  solve  the  equations 

da  db  dc 

Since  these  equations  could  not  be  solved,  estimates  o',  6',  c'  of  o,  6,  c  were  obtained 
and  the  function  expanded  by  a  Taylor  series  about  this  point,  namely 

/(o'  4-  a,  5'  +  /?,  c'  -i  y ,  x)  =  /(o',  6':  c',  x)  +  ^  0^  +  7  ^]/(® 

plu*  higher  order  terms  which  may  be  neglected  if  a,  ^  and  y  are  smaU ;  that  is.  if  the 
estimates  o',  6'  r.nd  c'  good.  The  normal  equations  are  then  derived  as  follows: 
Lot 

hence  ^ 

V,  "  ^  +  ^  •  0  for  a  minimum, 

l»t 

and  similarly  lor  dFIdb  and  dFfde. 

Applying  tite  approximate  value  from  the  Taylor  expansion,  we  have 


«.  etc. 
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which  gives  the  normal  equations 


a£| 

a£| 

These  equations  are  solved  to  give  values  of  a,  y  which  are  added  to  the  estimates 
of  n,  b  and  c  to  give  new  estimates,  and  the  process  may  l)e  iterated  until  a  minimum 
value  of  the  standard  deviation  is  obtained  to  a  desired  degree  of  accuracy.  For  this 
iterative  proces.s  to  converge  it  is  necessary  for  a,  ft  an<i  y  to  be  small,  that  is, 
good  initial  estimates  of  the  parameters  must  be  olitaiiuHl.  This  it«*rative  technique 
was  used  to  lit  the  observed  data  to  the  equation 

F  =  !;:( 1  -  y?/!)  e-*' +  a  In  ( I  + /yt). 
and  thus  find  the  best  fit  values  of  1^,  a,  ft  and  a. 


6.  Results 

The  initial  aims  of  the  studies  de  scribed  here  were  to  develop  a  method  for  the 
measurement  of  the  air  velocity  in  blast  waves,  aid  subsequently  to  use  this  tech¬ 
nique  to  study  the  air  motion  in  blast  waves  from  a  vdde  variety  of  charge  sizes. 
It  had  been  intended  to  make  these  measurements  in  the  range  of  intermediate 
shook  strengths  from  about  7*5  to  0-5  atm  overpressure.  At  high  shock  strengths, 
however,  the  temperature  behind  the  shock  was  such  that  the  sulphuric  acid 
droplets,  which  formed  the  trails,  were  vaporized.  This  effect  occurred  at  shock 
strengths  of  about  5  atm  overpressure  for  which  the  temperature  behind  the  shock  is 
300  °C.  Luohinskii  (1956)  gives  the  boiling-point  of  sulphuric  acid  formed  from  a 
water  vapour -f  sulphur  trioxido  mixture  as  being  in  a  range  from  280  to  330  ^C, 
depending  on  the  percentage  of  water  present.  The  temperature  behind  the  shock 
at  this  pressure  level  increases  with  distance  from  the  shock  front  (Brode  1959), 
which  would  account  for  the  observation  that,  in  the  limiting  region  for  this  effect, 
the  trails  do  not  disapiiear  at  the  shock,  where  the  overpressure  is  greatest,  but  at 
some  distance  betiind  the  shook.  The  majority  of  air  velocity  measurements  have 
therefore  been  restricted  to  the  region  of  shock  strengths  from  6  to  0-5  atm  over¬ 
pressure.  In  current  work  titanium  tetrachloride  is  bi>ing  used  as  a  tracer,  and  satib- 
factoiy’  results  have  been  obtained  at  higher  pressures. 

The  air  velocity  w  as  measured  in  blast  waves  produced  by  surface  burst  charges 
of  t.n.t.  of  the  following  masses:  30,  60, 100«),  10000,  40000  and  2000001b.  All  the 
charges  were  made  from  cast  pure  t.n.t.  with  central  Tetr>’tol  boosters.  The  boosters 
accounted  for  about  0*2  of  the  charge  mass.  The  30  and  60  lb.  charges  were  cast 
as  solid  spheres,  v  htle  the  large  charges  were  built  in  the  form  of  hemispheres  from 
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12  in.  X  12  in.  x  4in.  cast  block8,eaohofabout32*6Ib.  The  manufactureofthecharges 
has  been  described  by  Pennie,  Philips  &  Holdsworth  (i960)  and  Philips,  Ditto  & 
Holdsworth  (i960).  In  all  experiments  using  block-built  charges,  probes  were 
placed  throughout  the  explosive  and  uniform  detonation  was  confirmed  in  all 
oases. 


Fioobk  6.  Comparisons  of  peak  particle  velocities  for  surface  hurst  t.n.t.  chargee  of  60  lb. 
(A);  10000 Ib.  (  +  );  400001b.  (O)  and  2000001b.  (•).  The  line  shows  the  particle 
velocity  calculated  from  the  shock  velocity. 

The  results  frcm  all  experiments  were  analyzed  by  the  methods  described  in 
Hs  and  4,  and  the  valueo  of  the  peak  velocity,  V,,  determined  by  the  best  fit 
equation  to  the  velocity-time  data  are  plotted  against  sealed  distance  in  figure  6. 
Also  shown  is  the  variation  of  peak  particle  velocity  determined  by  applying  the 
observed  shock  velocity  in  continuity  equations.  The  variation  of  velocity  with 
time  at  specific  distances  has  also  been  compared  with  Brode’s  calculation  and  a 
typical  comparison  was  shown  in  figure  6.  Several  significant  difl'erences  may  be 
observed  between  the  two  curves.  The  actual  rate  of  decay  of  the  velocity  is  slower 
than  that  indicated  by  Erode  and  the  arrival  of  the  secondary  shook  is  sooner,  but, 
due  to  the  extended  outward  fiow  from  the  charge,  it  occurs  in  the  positive  flow 


velocity  (Mach  units) 
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period.  There  is  an  indication  from  some  experiments  that  the  flow  velocity  may 
increase  slightly  prior  to  the  arrival  of  the  .  econdaty  shock,  but  this  has  not  been 
confirmed.  The  reason  for  this  difference  will  be  discussed  (see  §6).  The  secondary 
shock  occurs  in  the  blast  wave  from  a  chemical  explosion,  and  originates  as  an 
inward  facing  shock  produced  by  the  deceleration  of  the  contact  surface  between 
the  detonation  products  and  the  atmosphere.  This  shock  implodes  at  the  origin 
and  is  then  reflected  outward,  to  following  primary  shock,  as  shown  in  figure  8. 


PlovRB  7.  Comparison  of  the  time  variation  of  velocity  at  a  scaled  distance  of  5*8  (t./S 
from  siufaco  burst  t.n.t.  charges  of  601b.  (O):  400001b.  (•)  and  200  0001b.  (A).  The 
line  is  the  least  sqiuires  tit  to  the  200000  lb.  data. 


Comparisons  were  also  made  between  the  velocity-time  variations  obtained  using 
the  different  explosive  yields,  and  tyiHoal  comparisons  at  a  specific  scaled  distance 
are  shown  in  figure  7.  At  no  distances  could  any  significant  differences  be  observed 
and  the  data  for  all  the  experiments  could  therefore  lie  combined  and  fitted  to 

i;(l-/»)e  *'-holn{l+/»). 

The  resulting  best  fit  values  of  V^,  <x,  fi,  a  and  t*  are  given  in  table  1,  where  I*  is 
the  time  for  the  velocity  to  fall  to  P'/e.  This  latter  parameter  is  given  since  it  has 
mmre  physical  significance  than  1//?  which  should  not  be  considered  as  the  positive 
dt  ration  in  this  case  but  rather  as  a  normalising  factor  for  time. 

All  of  the  information  derived  by  the  smoke  trail  experiments  may  be  presented 
on  a  single  z-t  diagram,  as  shown  in  figure  8.  This  information  includM  Uio  distauce 
versus  time  of  the  primary  and  secondary  shocks,  particle  trajectoi  iM  as  traced 
by  individual  smoke  elements,  the  isotachs  joining  points  of  equal  velocity  and  the 
distance  variation  of  f*. 


376 


J.  M.  Dewey 


Tabls  1.  Best  nr  parametbrs  of  equation  V  -^t)+a  In  (l+fit) 


peak 


scaled 

velocity. 

distance 

V, 

(ft./1b.*) 

(Mach  units) 

at 

P 

a 

<• 

80 

102 

0-820 

1-242 

0-848 

0-83 

8-5 

0-88 

0-460 

1-174 

0-500 

0-89 

60 

0-78 

0-418 

1-162 

0-467 

0-97 

6-8 

0-69 

0-374 

1-140 

0-437 

1-05 

70 

0-61 

0-343 

1-131 

0-410 

1-12 

7-8 

0-88 

0-318 

1-118 

0-384 

1-18 

80 

0-80 

0-300 

1-101 

0-360 

1-23 

8-6 

0-46 

0-290 

1-070 

0-332 

1-27 

0-0 

0-42 

0-279 

1-005 

0-313 

1-32 

08 

0-39 

0-268 

1-043 

0-294 

1-39 

too 

0*36 

0-288 

1-030 

0-276 

1-44 

10-8 

0-33 

0-246 

1-020 

0-264 

1-81 

110 

0-31 

0-238 

1-007 

0-252 

1-56 

118 

0*29 

0-228 

0-998 

0-239 

1-60 

120 

0-27 

0-220 

0-989 

0-229 

1-63 

12-8 

0-26 

0-211 

0-980 

0-220 

1-66 

130 

0*24 

0-206 

0-975 

0-210 

1-68 

PiouMi  8.  r-4  dMRntti  ■howing  the  trajeetorica  of  the  primary  (•>— )  and  areomhuy  ( -) 

■hock*:  the  particle  trajoctorir*  ( - );  the  inotach*  (  -  |  joining  point*  of  equal  Macdi 

velocity  and  f*(0)t  the  time  at  which  the  velocity  ha*  dece^^  to  I  -t  of  it*  p^  value. 
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6.  Discussion 

The  technique  for  the  measurement  of  air  velocity  in  blast  waves  has  been 
developed  to  a  stage  at  which  it  can  be  applied  to  firee-held  measuiements  on  a 
routine  basis,  and  with  inherent  errors  not  larger  than  those  of  other  blast  measure¬ 
ments.  The  smoke  is  observed  inside  the  shock  sphere  so  that  light  from  the  trails 
is  refracted  continuously  in  its  passage  through  the  high-density  region  behind  the 
shook,  and  also,  in  the  opposite  direction,  at  the  front  itself.  This  refraction  will 
give  a  slight  shift  to  the  image  of  the  trail,  but  calculations  based  on  an  assumed 
density  variation  behind  the  shock  indicated  that  the  apparent  displacement  would 
be  less  than  0*  1  %  in  the  region  of  interest  and ,  in  fact,  only  0*3  %  aii  a  shock  strength 
of  200  atm.  This  point  was  also  checked  experimentally  by  attempting  to  detect 
any  apparent  shift  of  the  marker  boards  as  they  were  *;aveloped  by  the  wave. 
No  such  a  shift  could  be  measured. 

The  peak  velocities  and  the  variation  of  velocity  with  time  at  specific  distances 
from  the  charge  have  been  compared  with  the  predictions  for  a  t.n.t.  explosion, 
given  by  Erode  (1959),  who  based  his  calculations  on  an  energy  release  of  1  kcal. 
He  assumed  a  figure  of  252kcal/mole  for  t.n.t.  with  a  density  of  l-Sg/cm*,  which 
gives  a  scaling  factor  of  9*4  fm*  a  1  lb.  chaige.  A  random  selection  of  813  cast  t.n.t. 
blocks  from  those  used  to  build  the  charges  for  the  experiments  described  above, 
had  a  mean  density  of  1*57  g/cm*.  For  cast  t.n.t.  of  this  density,  CotA  (1958)  quotes 
an  energy  yield  of  l*16kcal/g  which  gives  a  scaling  factmr  of  0*6  applicable  to 
Brode’s  results  for  a  1  lb.  charge.  Measurements  of  peak  (Nessure  have  been  made 
by  Groves  (personal  communication)  fw  a  large  number  of  air-burst  spherical 
t.n.t.  charges,  using  a  photogrammetrical  method  of  determining  the  shock 
velocity.  The  (uessure  distance  curve  from  these  experiments  is  shown  in  figure  9. 
A  comparison  wit'i  corresponding  results  by  Erode  gives  a  scaling  fMtor  of  10*0,  and 
this  factor  has  therefore  been  used  for  ail  the  air  velocity  oomparisoiu.  It  would 
appear  that  the  energy  >ield  of  t.n.t.  detonated  in  free  air  may  be  about  4  %  greatMr 
than  that  indicated  by  Cook. 

The  air  velocity  measurements  described  here  a-ere  made  with  surface  burst 
charges,  and  to  compMc  these  results  with  prediotkms  for  a  free  air  spherical 
explosion,  it  was  necessary  to  apply  a  ground  reflexion  factor.  This  was  done  by 
comparing  the  shock  overpressure-distance  relationships  of  air  and  ground  burst 
charges  and  finding  the  ratio  of  the  distances  at  which  a  specific  pressure  was  ob¬ 
served.  ITie  culie  of  thb  ratio  gives  a  reflexion  factor  which,  for  a  perfectly  reflecting 
surface  would  be  two.  Figure  0  shows  such  a  comparison  which  yields  reflexion 
factors  that  vary  in  the  region  of  present  int<mst  from  1*5  to  1*76,  depending  on 
the  distance  from  the  charfp'.  The  comparisons  between  the  observed  and  predicted 
velocities  have  been  made  using  the  appropriate  reflexion  factor  calculated  in  this 
way.  To  confirm  the  validity  of  this  prooedorB  the  air  velocity  sras  measured  in  the 
blast  waves  from  spherical  t.n.t.  60  lb.  charges  detonated  at  45  ft.  above  the  ground. 
Comparisons  between  the  free  air  and  surCMS  bursts  results  sealed  in  this  way,  are 
riiown  in  figure  10  and  no  significant  difference  can  be  observed. 

Whmi  it  became  clear  that  the  decay  of  air  velocity  in  Mast  waves  from  t.n.t. 
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expIoHions  wa8  not  following  the  expected  form,  comparisons  were  made  between 
the  observed  velocity  decay  and  that  derived  from  Brode’s  calculations,  f  Such  a 
comparison  was  shown  in  figure  6.  The  difference  between  the  observed  and  pre¬ 
dicted  velocities  was  studied  and  found  to  yield  a  consistent  relationship  of  the  form 
Av  =:  aln  (1  +fit)  where  Ar  is  the  velocity  difference  (see  figure  11).  Combination 
with  the  modified  Friedlander  equation  led  directly  to  the  form  quoted  above.  It 
must  be  pointed  out  that  this  equation  is  empirical  and  no  physical  significance 
can  be  attached  to  any  of  the  fitted  parameters,  other  than  V„  the  peak  velocity. 


Floras  •.  Variation  of  peak  o%-erpreMun«  with  dwtanoe  from  the  charge.  Solid  line  is  the 
free  air  relatiotiship.  The  broken  line  m  a  relatioiiahip  for  a  surface  burst  (dtarge.  The 
numbers  against  the  arrows  are  reflexion  factors  at  the  specified  pressures.  The  points 
(O)  are  derived  from  Brode’s  calculations  for  a  scaling  fsctor  of  10-0. 

It  la  suggested  that  the  departure  of  the  velooity  decay  from  that  of  a  olaseica) 
blast  wave  is  due  to  'after-burning'  iff  t-.n.t.  Trinitrotoluene  has  s  chemical  com¬ 
position  of  the  form  CyH»N,0,.  Cook  (1958)  gives  the  detemaiion  producta  of  1*59 
density  t.n.t.  ss  carbon  monoxide  (I).  oarl^  dioxide  (10),  water  (1),  nitrogen  (6-0), 
ammonia  (0'4),  methyl  aloohd  (4-5).  hydrogen  cyanide  (1*1)  and  oarbem  (14*2) 
wbne  the  figuree  in  brackets  are  the  number  iff  moloa  per  kilogramme.  Although  Uie 
exact  composition  of  these  producta  will  depend  on  the  nature  of  omifinement  of  the 
expkMve  (Jones  A  Miller  1948),  a  simplified  analysis,  assuming  stoichkNnetric 
mixtures  indicates  an  oxygen  deficiency  of  greater  than  60  %  by  mam.  On  several  of 
the  large-scale  experimenta  radiation  measuremonta  wore  made  to  determine 
fireball  black-body  tenperatmres,  such  as  thosedescribedbyTate  APiattman(i96a) 

t  The  ivpoii  of  Brodc's  calculations  for  t.n.t.  quotes  the  variation  of  bIset  psrameten  with 
distsneeftw  the  explosion.  St  specific  tirass.  For  prscticsl  studies  of  blsatwnvss  it  is  usually 
more  con^'cnient  to  know  the  variation  with  time  at  spceiflc  loestioM.  Dr  Brode  wm  kind 
enough  to  icml  hi*  original  computer  outputs  to  tho  author,  so  that  the  required  infnrmation 
could  be  derived  ami  used  in  the  oorapaiisons  disemued  here. 
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for  a  100 ton  t.n.t.  explosion.  These  measurements  showed  a  peak  temperature  of 
8600  °K  at  0-7  ms  after  detonation,  which  may  be  associated  with  the  shock  front 
temperature.  As  the  shock  front  expanded  and  cooled  it  became  transparent  to  the 


Kioi'rr  10.  C'ompariMin  of  tho  velooity>tinie  vsriatioiu  from  froe  air  and  aurfaoe  butat  60  lb. 
t.n.t.  ohargea.  + ,  Free  air  roeaauremeni  at  a  aoaled  diatanoe  of  4*8  h./Si  0>  surface 
burst  ineasurmient  at  a  scaled  distance  of  5-6  line  shows  the  least  squares  fit 
to  tho  eoinbinod  data. 


FmtrM  II.  Difihfenee  betsreo  thn  pmlictad  (Brode)  and  obssnrwl  partids  veloeitiea 
pbUed  againsl  In  (I  -f  t'r).  r  is  Uw  dumtion  of  tha  [iwdiotsd  posWve  flow. 
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radiation  from  the  detonation  products,  yielding  a  black- body  temperature  of 
about  3000  °K.  These  products  continued  to  radiate  for  several  seconds  after  detona¬ 
tion  at  a  temperature  of  about  2500  °K.  It  iherefore  appears  that  the  detonation 
products  are  of  a  composition  and  temperature  such  that  several  exothermic  reac¬ 
tions  will  take  place  in  the  presence  of  atmospheric  oxygen,  thus  serving  as  a 
relatively  slow  but  significant  enei^  source  for  some  time  after  detonation.  It  is 
suggested  that  this  additional  release  of  energy  accounts  for  the  non-classical  decay 
of  velocity  within  the  blast  wave. 


Kko'sk  lj.  C^npsriwm  of  the  obsnn'od  mmI  pmltrtMi  contsH  surisco  tfajoptorini  for  a 
XOOOOCI  lb.  t.n.t.  aurfare  iMiTHt  clMr{p\  .Solid  lino  ahowa  tho  ooaen'od  trajeotory  and  tho 
brokm  Urn*  the  trajootorx*  po’dirtod  from  Brodo’a  ralmilationa.  Points  aro  the  obacwod 
(•)  ami  pmlictod  main  ahopk  trajortorios. 

Thia  effect  is  further  illustrated  by  a  comparison  between  tin*  olwerved  growth 
of  the  fireball  and  the  trajectory  of  the  contact  surface  predicted  by  firode.  8ueh 
a  comparison  is  shown  in  figure  12.  It  arill  be  seen  that  there  is  excellent  agreement 
between  the  two  shock  trajectories,  indicating  that  the  correct  scaling  factor  has 
been  used,  but  a  consklerable  difference  between  tlie  trajectories  of  the  two  contact 
surfaces  indicates  a  process  such  as  'after-burning*  for  which  no  idlowance  was 
made  in  Brode's  calculations.  This  phenommwn  also  has  been  observed  when 
exidoaivee  of  ktw  oxygen  balance  are  used  for  shallow  blasting.  When  the  com- 
biMtkm  ]«rodacts  escape  and  mix  with  the  atmosphere,  a  secondary  explosion 
ocean. 

In  order  to  test  this  *  after -burning  *  hypothesn,  it  would  be  deeirable  to  study  the 
Mast  waves  produced  by  t.n.t.  charges  detonated  in  an  oxjgeo-free  aUnoaphere. 
However,  pore  t.n.t.  carmot  be  detonated  easily  in  msssps  Imb  than  thtm%  8lb., 


(Ifaeh  units) 
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and  facilities  for  performing  suoh  an  experiment  ve  not  at  present  available.  An 
alternative  was  to  make  measurements  in  blast  waves  produced  by  explosives  with 
a  high  oxygen  balance,  such  as  nitroglycerine.  Cook  (1958)  gives  the  detonation 
products  of  nitroglycerine,  with  a  specific  gravity  of  1*6,  as  carbon  dioxide  (13*2), 
water  (11),  nitrogen  (S*9),  nitric  oxide  (1*4)  and  oxygen  (0*4)  where  the  figures  in 
brackets  indicate  moles  per  kilogramme.  The  high  sensitivity  of  this  explosive  makes 
it  difficult  to  use  for  experimental  purposes,  but  a  series  of  experiments  was  carried 
out  with  601b.  spherical  charges  of  Qeogel,  an  explosive  produced  by  Canaciian 


Fiuuaa  IS.  Oomparisoa  of  the  velovity-time  variatiaas  at  a  MaM  dwtaaoe  nf 
5*8  ti.fS  flraia  80  lb.  aurfiMW  bunt  ebaifas  of  and  Ueofel. 

Industries  Ltd.  and  composed  of  50%  nitroglyoerhte.  'Hie  manuhtoturm  quoted 
an  MMfgy  yield  1  %  hi^ier  than  tn.t.  and  a  3*5%  excess  trf'oxygm  in  the  explosive 
products.  The  prearars*distanee  curves  obtained  for  air  bunt  and  ground  burst 
CSeogel  charges  had  a  different  slope  from  thnae  oX  t.n.t.  urith  a  oroas-over  point  at 
about  I'TSaUn.  oonreaponding  to  a  distance  about  6ft.  from  a  1  lb.  charge. 
Closer  to  the  charge  the  pressures  from  t.n.t.  were  greater  and  at  longer  distancas 
they  wtwe  ie«  titan  those  for  Geogel.  This  was  to  have  been  expected  from  the 
slower  detonation  veloelty  of  Geogel  quoted  at  6300m/s  compart  urith  6e0Om/s 
for  t.n.t.  For  convenience,  comparisons  betweiitn  the  air  velooities  in  the  two  blast 
waves  were  made  riose  to  this  cross-over  point.  F%nn  13  shows  a  comparison  be¬ 
tween  the  veloeity-time  variations  far  t.n.t.  and  Oeogel.  Ihere  k  no  measurable 
difference  between  the  two  decays  until  the  vriocHy  has  fitOen  to  about  Mach  0*t, 
from  which  point  the  Qeogel  curve  continues  to  foO  in  the  daseieal  maonsr,  givinf 
a  sealed  positive  duration  of  about  2*3.  The  difference  between  the  blast  waves  from 
the  two  explosives  is  shown  even  mniv  dearly  by  the  comparison  between  the 
particle  trajeetories  shown  in  figure  14. 
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In  a  series  of  ezperime'its  with  601b.  and  Geogel  charges,  Uie  detonationa 
were  photographed  at  64  frames/s  on  a  colour  film.  A  luminous  fireball  £nnn  the 
t.n.t.  exploskms  existed  for  at  least  600  ms,  but  for  less  than  17  ms  few  the  Geogel. 

It  was  recently  brought  to  the  author’s  attention  that  a  limited  number  of  smoke 
puffs  and  trails  had  been  used  on  some  of  the  early  nuclear  explosions  detonated  by 
the  United  States,  and,  on  request,  a  copy  of  a  100  frames/s  film  showing  the  smoke 
movement,  was  supplied  by  the  U.S.  Government.  The  small  number  of  smoke 
trails  that  had  been  used  did  not  make  possible  a  complete  particle  velocity  analysis. 


14.  OmqMinMiiw  of  tlw  {wirtk4»  in  btasl  bom  M  lb.  muCsco 

bum  vbariito  of  ( o )  awi  ( -f ). 

but  sewtal  partide  trajectories  oaukl  be  followed.  Analysis  of  the  shock  velocily 
from  the  nndear  explooloa  indicated  a  blast  yield  of  dxwt  63  kilotons.  This  figure 
was  used  to  scale  the  results  so  that  they  could  be  compared  with  t.n.t.  d^a.  and 
Mich  a  comparison  b  shown  in  figure  16.  Tbs  initial  scaled  particle  dnplaoements 
from  the  tiro  sxploaioaB  are  very  Mmilar,  but  differ  signifioantly  at  later  times  when 
the  nuclear  Mast  deoays  in  the  classical  fonn.  ilddi  curves  wars  derived  from  free 
air  eqiloaiona  and  diecontinuitieo  are  caused  by  the  arrival  of  the  shodt  reSected 
from  the  ground. 

Taylor  (1950)  pointed  out  that  a  major  difikidty  fan  the  analyris  of  i^ph«rioal 
blasts  is  that  the  shook  wave  leaves  the  sir  in  a  state  in  whkh  the  entropy  deoreases 
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radially,  so  that  after  its  passage,  when  the  air  has  returned  to  atmospheric  pressure, 
the  air  temperature  decreases  with  increasing  distance  from  the  centre  of  the  ex¬ 
plosion.  For  this  reason  the  density  is  not  a  single  valued  function  of  the  pressure 
in  a  blast  wave,  and  if  the  pressure-time  relation  is  measured  at  a  point,  it  is  not 
possible  to  derive  directly  the  density -time  relation  unless  the  particle  trajectories 
are  known.  In  this  case  the  trajectory  of  the  air  which  was  monitored  by  the  pres¬ 
sure  transducer  at  uiy  instant  can  be  traced  back  to  its  position  when  it  was 
shocked.  The  Rankine-Hugoniot  equation  can  then  be  applied  at  that  point  and  the 
density  at  the  original  position  and  time  calculated  from  the  adiabatic  equation. 


Fiovas  15.  Contparwow  the  particle  trajectorim  in  blart  wa\-e«  from 
air  bunt  t.n.t.  (£i)  and  aurli«r  (#)  rxfikMiowi. 


Sttdi  a  oonverskm  has,  in  fact,  been  done  with  reasonable  success  (Dewey  ft  Anson 
1963).  On  the  other  band,  if  the  particle  trajectories  are  known,  it  is  possible  to 
oakuiate  the  density  directly  in  terms  of  the  ambient  value  by  applying  the  con¬ 
servation  of  a  mam  equation  w  hich  is  unaffected  by  the  entropy  change  at  the  shock. 
The  equation  may  be  written  for  a  sphcrkal  flow  in  the  form 

p  ^ 

^  “  r*  {9r!?xy 

where  is  the  aml^t  density  and  p  the  density  of  an  air  element  at  a  tiUH'  f 
and  radidl  poeition  r,  the  original  poeition  of  whk^  was  r.  The  partial  derivation 
dr/dbrmay  be  detormioed  from  the  observed  relationship  between  r  and  X  at  aoonstaat 
time.  A  family  of  sod)  curves  for  a  SOOOOOtb.t.n.t  exploeion  is  shown  in  figure  16. 
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The  hydrostatic  pressure  may  now  be  determined  from  the  density  by  using  the 
adiabatic  equation  as  described  above.  The  velocity  of  sound  and  thus  temperature 
ratios  can  then  be  found.  This  would  be  of  great  advantage  since  it  would  enable 
velocity  methods  to  be  applied  to  redected  and  secondary  shocks.  The  application  of 
the  above  relations  requires  the  assumption  of  a  constant  ratio  of  specific  heats 
throughout  the  wave,  and  this  is  probably  valid  for  blast  waves  of  strengths  less 
than  about  seven  atmospheres  overpressure.  These  relations  have  been  applied  in 


Fiocrb  16.  The  radial  positions  {r)  of  elements  in  a  blast  wave  plotted  against  their  original 
positiors  (f)  at  the  stated  times  (t).  Data  obtained  with  a  200000  lb.  surface  burst  charge. 
The  slopes  ol  these  ciu^es  arc  a  rneasiwe  of  the  radial  compression  from  which  the  density 
ratio  can  bo  calculated. 

preliminary  hand  calculations  which  showed  good  agreement  with  the  pressures 
and  densities  measured  independently.  Such  an  analysis  is  well  suited  to  a  digital 
computer  and  would  permit  a  complete  description  of  all  the  physical  parameters 
within  a  blast  wave  in  the  range  discussed,  using  only  photographic  instrumentation. 

7.  Conclusions 

The  smoke  trail  displacement  technique  has  enabled  the  particle  velocities  to  be 
measured  in  blast  waves  from  t.n.t.  explosions  varying  in  mass  from  30  to  200  000  lb. 
Throughout  this  range  the  results  have  scaled  completely,  using  a  factor 

s^nwpjp), 

where  If'  is  the  charge  mass,  P  the  ambient  pressure  and  a  standard  pressure. 
Good  scaling  was  also  achieved  between  air  burst  and  ground  burst  explosions  by 
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using  a  ground  reflexion  factor  in  the  range  of  1-5  to  1*76,  depending  on  the  distance 
from  the  charge.  There  was  no  significant  difference  between  the  measured  velo¬ 
cities  close  to  the  shock  and  those  predicted  from  Brode’s  (1959)  calculations 
after  applying  a  scaling  factor  of  10-0  for  a  1  lb.  charge.  At  later  times  the  decay  of 
velocity  did  not  follow  the  classical  form  and  there  was  an  extended  outw'ard  flow 
from  the  charge.  The  shape  of  the  velocity  decay  could  be  described  by  an  empirical 
equation  of  the  form  ^  ^  f;e-'ii-^)+<.ln(I+^). 

It  is  postulated  that  the  extended  flow  is  due  to  the  after-burning  of  the  t.n.t. 
detonation  products  in  the  presence  of  atmospheric  oxygen,  and  it  has  been  shown 
that  this  phenomenon  does  not  occur  for  an  explosive  with  a  high  oxygen  balance 
or  for  a  nuclear  explosion.  It  has  been  shown  that  for  the  range  of  shock  strengths 
studied,  it  is  possible,  in  theory,  to  determine  the  time  and  spatial  variation  of  all 
the  physical  parameters  of  blast  waves  from  the  photographic  analysis  of  the  par¬ 
ticle  trajectories.  These  analyses  will  be  the  object  of  future  studies. 
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Field  Section  under  the  leadership  of  Mr  N.  Spackman,  the  Munitions  Section  under 
Mr  W.  Ditto  for  provision  of  the  explosive  charges,  black  smoke  background  and  the 
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